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ABSTRACT 

We report extensive high-resolution spectroscopic observations and U-band differential photometry 
of the slightly eccentric 7.02-day detached eclipsing binary VSOlMon (A6m-|-F0), which we use to 
determine its absolute dimensions to high precision (0.3% for the masses and 1.8% for the radii, 
or better). The absolute masses, radii, and temperatures are Ma = 1.6455 ± 0.0043 Mq, Ra = 
1.888 ± 0.029 Rq, and = 7510 ± 100 K for the primary and Mb = 1.4588 ± 0.0025 Mq, R^ = 

1.592 ± 0.028 Rq, and = 7000 ± 90 K for the secondary. Apsidal motion has been detected, to 

which General Relativity contributes approximately 70%. The primary star is found to be a metallic- 
line A star. A detailed chemical analysis of the disentangled spectra yields abundances for more 
than a dozen elements in each star. Based on the secondary, the system metallicity is near solar: 
[Fe/H] = -1-0.01 ± 0.06. Lithium is detected in the secondary but not in the primary. A comparison 
with current stellar evolution models shows a good match to the measured properties at an age of 
about 1.1 Gyr. 

Subject headings: binaries: eclipsing — stars: evolution — stars: fundamental parameters — stars: 
individual (V501 Mon) — techniques: photometric 


1. INTRODUCTION 

The eclipsing binary V501Mon (with alternate names 
2MASS J06404172-0106400 and TYG 4799-1943-1; V = 
12. 32, SpT A6m-|-F0 as determined here) was discovered 
by IWachmannI (|1966f) in a photographic survey, and was 
given the original designation HBV 444. He proposed 
an orbital period of 7.021115 days, near the value we 
find below, and found a displaced secondary minimum 
indicating an eccentric orbi t. A reexamination of th e 
original plates was done bv iBossen fc Klawitted (1197211 . 
essentially confirming the original results, but with bet¬ 
ter precision. Since then, numerous times of eclipse have 
been published, and V501 Mon has been mentioned as 
a target worthy of followup for apsidal motion studies, 
possibly featuring an important c ontribution to the pre ¬ 
cession from General Relativity (|GimenedI1985L [l995ll . 
[Sebastian et al.l (|2012ll classified V501Mon as of spec¬ 
tral type A9 III. To date there has been no spectroscopic 
study of the system. 

In this paper we present new differential photometry 
of V501Mon as well as extensive high-resolution spec¬ 
troscopy, the combination of which enables us to derive 
accurate absolute dimensions for the components as well 
as their chemical abundances for more than a dozen el¬ 
ements. As we describe below, the primary is found to 
be a metallic-line A star. Most such objects are mem- 
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bers of binary systems (|Abtl I196II : lAbt fc Levvl 1198511 . 
Among the dozen or so Am binaries with well deter¬ 
mined properties, only a handful have been subjected to 
detailed abundance analyses fseelLvubimkov et al.llI996l : 

[Torres et al.ll20I2l : [Pavlovski et al.[[2014ll : thus V50I Mon 

adds significantly to this class of objects, and may assist 
in understanding the Am phenomenon. 

2. ECLIPSE TIMINGS 

Measurements of the times of eclipse for V50I Mon col¬ 
lected from the literature are listed in Table [U and cover 
approximately seven decade s. An ephe meris curve so¬ 
lution using the method of [LacH ([I99211 gave clear in¬ 
dicati ons of apsidal motion, as suspected by [Gimenej 
([1985I1 . and resulted in a best fit eccentricity of e = 
0.I3I3 ± 0.0019, longitude of periastron w = 231?54 ± 
0?71, and reference epoch of primary minimum HJD 
2,453,401.65029 ± 0.00060, where the inclination angle 
has been held fixed at the value determined from our 
light-curve fit below (i = 88?022). The apsidal period, 
however, is still very uncertain (U = 16,000 ± 11,000 yr), 
and strong correlations are present between the eccentric¬ 
ity and the rate of apsidal motion in this “unconstrained” 
solution. Given that the spectroscopic measurements we 
present below provide a stronger handle on e, a preferable 
approach is to perform a simultaneous fit of the eclipse 
timings and the radial velocities, which should constrain 
the parameters better. We defer a discussion of this so¬ 
lution until Section m 

3. SPECTROSCOPIG OBSERVATIONS AND RADIAL 
VELOCITIES 

Spectroscopic observations of V501 Mon were carried 
out with three different instruments. They began at the 
Harvard-Smithsonian Genter for Astrophysics (CfA) in 
2005 November, u sing the now d ecommissioned Digital 
Speedometer (DS; [Lathaml[l992[) mounted on the 1.5 m 
Tillinghast reflector at the Fred L. Whipple Observatory 
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TABLE 1 

Times of eclipse for VSOIMon. 


HJD 

(2,400,000-1-) 

e 

(days) 

Eclipse 

Type 

(O-C) 

(days) 

Source 

30077.2220 


1 

pg 

-1-0.01762 

1 

30431.3800 


2 

pg 

-0.01654 

1 

31077.3560 


2 

pg 

-1-0.00819 

1 

31144.4230 


1 

pg 

-0.00467 

1 

31144.4280 


1 

pg 

-1-0.00033 

1 

31786.4970 


2 

pg 

-1-0.00703 

1 

31846.5810 


1 

pg 

-1-0.03275 

1 

33703.3080 


2 

pg 

-1-0.02785 

1 

33710.3090 


2 

pg 

-1-0.00764 

1 

34040.2580 


2 

pg 

-0.04021 

1 

34068.3530 


2 

pg 

-0.03005 

1 

34086.3060 


1 

pg 

-0.00692 

1 

34086.3400 


1 

pg 

-1-0.02708 

1 

34416.3070 


1 

pg 

-0.00259 

1 

34444.4320 


1 

pg 

-1-0.03758 

1 

34451.3920 


1 

pg 

-0.02362 

1 

34767.3870 


1 

pg 

-1-0.01711 

1 

34767.4250 


1 

pg 

-1-0.05511 

1 

34770.4820 


2 

pg 

-0.02199 

1 

34781.3920 


1 

pg 

-0.02030 

1 

34781.4220 


1 

pg 

-1-0.00970 

1 

34795.4210 


1 

pg 

-0.03371 

1 

35107.5130 


2 

pg 

-0.00905 

1 

35160.5080 


1 

pg 

-0.04941 

1 

35160.5390 


1 

pg 

-0.01841 

1 

52320.38766 

0.0025 

1 

ccd 

-1-0.00322 

2 

52983.5215 

0.0002 

2 

ccd 

-1-0.00009 

3 

53376.7091 

0.0005 

2 

ccd 

-0.00004 

4 

53390.7530 

0.0006 

2 

ccd 

-1-0.00144 

4 

53401.6502 

0.0005 

1 

ccd 

-1-0.00006 

4 

53404.7933 

0.0005 

2 

ccd 

-0.00068 

4 

53671.6010 

0.0005 

2 

ccd 

-1-0.00106 

5 

54422.8674 

0.0005 

2 

ccd 

-0.00195 

6 

54450.9531 

0.0007 

2 

ccd 

-0.00109 

6 

54454.8302 

0.0004 

1 

ccd 

-0.00081 

6 

54475.8934 

0.0008 

1 

ccd 

-0.00123 

6 

54787.9741 

0.0006 

2 

ccd 

-1-0.00185 

6 

54791.8490 

0.0010 

1 

ccd 

-1-0.00011 

6 

55135.8866 

0.0019 

1 

ccd 

-0.00137 

7 

55170.9960 

0.0013 

1 

ccd 

-1-0.00200 

7 

55493.9715 

0.0009 

1 

ccd 

-1-0.00203 

7 

55500.9939 

0.0011 

1 

ccd 

-1-0.00323 

7 

56339.65867 

0.00026 

2 

ccd 

-0.00088 

8 

56339.65954 

0.00038 

2 

ccd 

-0.00001 

8 

56339.66015 

0.00022 

2 

ccd 

-1-0.00060 

8 

56676.6771 

0.0009 

2 

ccd 

-0.00050 

9 


Note. — Measurement errors (e) are listed as published. 
“Eclipse” is 1 for a primary eclipse, 2 for a secondary eclipse. 
“Type” is pg for photographic measurements and ccd for more 
modern determinations. Uncertainties for the photographic mea¬ 
surements are estimated to be e = 0.025 days, as explained in 
Section |4] O — C residuals are computed from the combined 
fit described there. Sources for the times of eclipse are: (1) 
http://var2.ast ro.cz/EN/brno/eclipsing_binaries.php, (2) 

IBrat et al.l Il2()07l l: (3j Measured from a mi nimum observed b y 
M. Wolf ( priv. com m.); (4 ) |Lacyl 11200(3): (5) Zeida et al.l 1(200^^ 
(6) ILacvl 1120091) : (7) ILac^ 1120111) : (8) IPiethe ml 112013 ): (9) ILacvl 
1120141) . 


on Mount Hopkins (AZ). Seven spectra were recorded 
through 2009 March with an intensified photon-counting 
Reticon detector, and cover a narrow span of 45 A cen¬ 
tered at 5190 A (Mg Ib triplet). The resolving power of 
this instrument was R k, 35,000, and the signal-to-noise 
ratios of the spectra range from 13 to 22 per resolution 
element of 8.5 kms“^. Dusk and dawn exposures of the 
sky were taken nightly to monitor the velocity zero point. 

Thirty seven additional spectra were gathered from 
2009 November to 2015 February with the Tillinghast 


Reflector Echelle Spectrograph (TRES; iFiireszl[200811 on 
the same telescope. This bench-mounted, fiber-fed in¬ 
strument provides a resolving power of i? sa 44,000 in 
51 orders over the wavelength span 3900-9100A. The 
signal-to-noise ratios of the 37 spectra range from 8 to 56 
per resolution element of 6.8 kms“^. lAU radial-velocity 
standard stars were observed each night. 

Between 2011 October and 2015 February we also ob¬ 
tained 57 usable spectra of V501 Mon with the Ten¬ 
nessee State University 2 m Automatic Spectroscopic 
Telescope (AST) and a fi ber-fed echelle spectrograph 
([Eaton fc Williamsonll200^ at Fairborn Observatory in 
southeast Arizona. The detector for these observa¬ 
tions was a Fairchild 486 CCD, with 15 /rm pixels in a 
4096 X 4096 format. The spectrograms have 48 orders 
ranging from 3800-8260 A. Because of the faintness of 
V501Mon (V = 12.32), we used a fiber that produced a 
spectral resolution of 0.4 A, co rresponding to a r esolving 
power of 15,000 at 6000 A. See lFekel et al.l (|2013[ 1 for ad¬ 
ditional information about the AST facility. Our spectra 
have typical signal-to-noise ratios per resolution element 
of 40 at 6000 A. 

Radial velocities (RVs) from the CfA spectra were mea¬ 
sured using the t wo-dimensional cross-c orrelation tech¬ 
nique TODCOR (jZucker fc Mazet3ll994f) . with synthetic 
templates taken from a large library of calculated spec- 
tra based on model atmospheres by R. L. Kurucz (see 
[Nordstrom et al.l 11994 iLatham et al.l 120021) , restricted 
to the Mg Ib region. Template parameters (effective 
temperature, surface gravity, metallicity, and rotational 
broadening) were selected b y running gr i ds of cross¬ 
correlations as described bv [Torres et alJ (I2002D . We 
adopted initial surface gravities of log 5 = TO, close to 
our final values in Section |3 and solar metallicity. Fur¬ 
ther experiments showed that while solar composition is 
favored for the secondary (star B), a better match to the 
observed spectra is achieved by adopting a metallicity for 
the primary (star A) of [Fe/H] = -1-0.5, which we have 
often seen in the past in objects with chemical peculiari¬ 
ties typical of Am stars. This was a first indication that 
the primary component of V501 Mon may be metallic- 
lined. The temperatures for the primary and secondary 
that maximize the cross-correlation coefficient averaged 
over all exposures are 7280 K and 6950 K, although these 
are sensitive to the adopted metallicity and surface grav¬ 
ity and could be biased if the primary composition is 
anomalous. Formal uncertainties for these temperatures 
are 200 K each. The rotational broadenings were deter¬ 
mined to be usinf = 17±1 kms“^ and 11±2 kms“^, re¬ 
spectively. The heliocentric velocities we obtained from 
the DS and TRES spectra were placed on a common 
frame and are listed in Table [2] Those from the DS 
include corrections for potential sy stematic effects due 
to the narrow spectral window fsee ILatham et aI]ll996D 
determined from n umerical simulations as described by 
[Torres et al.l (|1997f) . These corrections are smaller than 
0.5 kms“^ for both stars, and have little effect on the 
results presented later. 

We determined the light ratio separa tely for the seven 
DS and the 37 TRES spectra following I Zucker fc Mazehl 
(|1994f) . and obtained values of £b/-^a = 0.58 ± 0.06 and 
0.53±0.02, respectively, at a mean wavelength of 5190 A. 
We point out, however, that these values are also sen- 
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TABLE 2 

Heliocentric radial velocity measurements of V501 Mon from CfA. 


HJD 

(2,400,000-1-) 

RVa 

(kms“^) 

eA 

(kms“^) 

(O - C)a 
( kms“l) 

RVb 

(kms~^) 

eB 

(kms“^) 

(O - C)b 
( km s“^) 

Orbital 

phase 

Instrument 

53690.9255 

-84.96 

0.64 

-0.42 

84.73 

2.32 

+3.61 

0.2002 

DS 

54048.9456 

-82.16 

0.47 

+0.50 

76.84 

1.72 

-2.16 

0.1915 

DS 

54070.9780 

-76.65 

0.39 

+0.10 

72.80 

1.40 

+0.47 

0.3295 

DS 

54136.8728 

62.18 

0.46 

-0.45 

-85.06 

1.66 

-0.16 

0.7146 

DS 

54424.9270 

60.25 

0.37 

-0.07 

-82.04 

1.32 

+0.25 

0.7409 

DS 


Note. — This table is available in its entirety in machine-readable and Virtual Observatory (VO) forms in 
the online journal. A portion is shown here for guidance regarding its form and content. 


sitive to the metallicity and temperatures adopted for 
the cross-correlation templates, and may be biased if the 
composition of the primary is peculiar (which could in 
turn affect the temperature derived for this star). For 
example, using a primary temperature of 7500 K, closer 
to the value found below, reduces the light ratio from 
TRES to -^b/^a = 0.49 ± 0.02. Although a template 
mismatch may affect the derived temperatures and light 
ratio, experience has shown that it has little effect on 
V sinf or on t h e radi al velocities. 

iFekel et'all (I2009D gave a general explanation of the 
velocity measurement from our Fairborn echelle spec¬ 
tra. For V501 Mon we have used our solar-type line list 
and measured velocities for the lines in that list that 
are just in the orders that cover the wavelength region 
4920-7100 A. Our velocities were determined by fitting 
the individual lines w ith rotational broadening functions 
(jLacv fc Fekell 1201 Ih that allowed both the depth and 
width of the line fits to vary. Our unpublished measure¬ 
ments of several lAU solar-type velocity standards show 
that these Fairborn Observatory velocities have a zero- 
p oint offset of —0.6 kms“^ when compared to the results 
of iScarfd (j2010D . Therefore, 0.6 kms“^ has been added 
to each velocity. We list the final values in Table [H 

From the solar-type line list for the AST spectra, the 
average projected rotational velocities of components A 
and B are 16.5 and 12.4 kms“^. These vsini values have 
an estimated uncertainty of 1 kms“^. 

4. SPECTROSCOPIC ORBITAL SOLUTION 

Separate orbital solutions using the CfA and Fairborn 
velocities are listed in Table IH and are generally con¬ 
sistent with each other. In order to optimally combine 
the information contained in the radial velocities and in 
the eclipse timings, we carried out a joint solution us¬ 
ing the two RV data sets and the times of eclipse from 
Section [51 solving for the spectroscopic orbital elements 
and the apsidal motion simultaneously. The pre scrip¬ 
tion f or the latter is based on the algorithm of iLacvl 
(jl992D . The long baseline of the eclipse timings helps 
to constrain the period, and the spectroscopic measure¬ 
ments constrain the orbital eccentricity, which is other¬ 
wise strongly correlated with the apsidal motion if using 
only the timings to determine doj/dt. Formal measure¬ 
ment uncertainties were used to establish weights, and 
adjusted iteratively during the solution so as to achieve 
reduced chi-squared values near unity, separately for the 
primary and secondary components and for each type of 
observation. For the photographic times of eclipse that 
have no published uncertainties our procedure assigned 
an average error of 0.025 days. Photoelectric/CCD tim- 



Phase 

Fig. 1.— Spectroscopic orbital solution for V501 Mon along with 
the measurements. Primary velocities are shown with filled sym¬ 
bols, and the dotted line represents the center-of-mass velocity. 
The bottom panels display the residuals, with those from Fairborn 
Observatory displaced vertically for clarity. 

ings with published errors in Table [T] were scaled by fac¬ 
tors of 1.7 for the primary and 2.2 for the secondary. 
The inclination angle was held fixed according to the 
value determined from our light curve fit in Section |5j 
Tests show that its small error contributes negligibly to 
the uncertainties in all other adjusted quantities. 

The results of our weighted least-squares fit are given 
in the last column of Table U] Though improved over the 
preliminary value listed in Section |2 the apsidal motion 
period U remains relatively uncertain at 15400±8300 yr. 
For each RV set we included as an adjustable parameter 
an offset between the zero points of the primary and sec¬ 
ondary velocities, which for the CfA measurements may 
arise because of a template mismatch, or for other rea¬ 
sons. We also allowed for a possible overall shift between 
the Fairborn and CfA velocities. All of these shifts turn 
out to be small, but are statistically significant. We illus¬ 
trate the solution for the spectroscopic elements in Fig¬ 
ure [I] which shows the measurements and their O — C 
residuals from the fit, separately for the Fairborn and 
CfA data sets. The residuals are listed in Tables [5] and |31 
Those from the eclipse timings and the best-fit ephemeris 
curve are displayed in Figure |2j 
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TABLE 3 

Heliocentric radial velocity measurements of V501 Mon from Fairborn Observatory. 


HJD 

(2,400,000-1-) 

RVa 

(km s“^) 

EA 

(km s“^) 

(O - C)a 
( kms“l) 

RVb 

(kms“^) 

EB 

(kms“^) 

(O - C)b 
( km s“^) 

Orbital 

phase 

55864.835 

46.10 

0.69 

-0.36 

-65.50 

1.26 

-1-0.61 

0.8207 

55928.777 

14.40 

0.69 

-1-0.05 

-31.70 

1.26 

-1.81 

0.9277 

55945.678 

-74.40 

0.69 

-0.03 

71.40 

1.26 

-1-1.21 

0.3348 

55949.894 

13.40 

0.69 

-1-1.76 

-25.70 

1.26 

-1-1.14 

0.9353 

55958.876 

-86.60 

0.69 

-1-0.03 

85.70 

1.26 

-1-1.68 

0.2145 


Note. — This table is available in its entirety in machine-readable and Virtual Observatory 
(VO) forms in the online journal. A portion is shown here for guidance regarding its form and 
content. 


TABLE 4 

Spectroscopic/ephemeris-curve solutions for V501MON. 


Parameter 

CfA 

Fairborn 

RVs + timings 

Asid (days). 



7.02120771 ± 0.00000097 

Panom (days). 

7.021201 ± 0.000014 

7.021254 ± 0.000061 

7.0212164 ± 0.0000053 

7 (kms~^). 

-6.91 ± 0.11 

-6.692 ± 0.099 

-7.067 ± 0.045 

i^A (kms“^). 

76.914 ± 0.054 

76.91 ± 0.12 

76.866 ± 0.047 

Kb (kms“^). 

86.58 ± 0.13 

87.22 ± 0.19 

86.71 ± 0.11 

e. 

0.1359 ± 0.0015 

0.1348 ± 0.0015 

0.13388 ± 0.00059 

UJA at Tjnin 1 (deg). 

232.19 ± 0.66 

233.07 ± 0.65 

232.43 ± 0.21 

Tlnin 1 (HJD-2,400,000). 

53401.6532 ± 0.0072 

53401.633 ± 0.025 

53401.65013 ± 0.00044 

Tperi (HJD-2,400,000). 

53404.224 ± 0.012 

53404.227 ± 0.057 

53404.2302 ± 0.0039 

dojldt (deg cycle" ^). 



0.00045 ± 0.00024 

u (yr) . 



15400 ± 8300 

* (deg). 



88.022 (fixed) 

Ma sin^ i (Mq) . 

1.6374 ± 0.0049 

1.6628 ± 0.0078 

1.6426 ± 0.0043 

Mb sin^ i {Mq) .. 

1.4545 ± 0.0029 

1.4663 ± 0.0054 

1.4562 ± 0.0025 

aA sin i (10® km). 

7.3570 ± 0.0045 

7.358 ± 0.011 

7.3545 ± 0.0043 

as sin 2 (10® km). 

8.282 ± 0.012 

8.344 ± 0.018 

8.296 ± 0.010 

a sin 2 (Rq) . 

22.481 ± 0.019 

22.571 ± 0.031 

22.498 ± 0.016 

q = Mb/Ma . 

0.8883 ± 0.0013 

0.8818 ± 0.0023 

0.8865 ± 0.0012 

ARV CfA (prim—sec) (kms“^). 

-0.50 ± 0.25 


-0.79 ± 0.13 

ARV Fairborn (prim—sec) (kms“^)... 


-0.60 ± 0.20 

-0.43 ± 0.18 

ARV (CfA-Fairborn) (kms-i). 



-0.426 ± 0.097 

C7A (kms~^). 

0.26 

0.64 

0.28 / 0.64 

C7B (km s“^). 

0.74 

1.08 

0.73 / 1.16 

Time span (days).. 

3372.9 

1207.9 

26995.5 

Kbv . 

44 

57 

44 / 57 

Aecl. 



46 


Note. — Symbols in the first column represent the sidereal and anomalistic periods, the center-of-mass velocity (which for the 
combined fit is on the Of A reference system), the velocity semi-amplitudes, eccentricity, longitude of periastron of the primary, 
reference time of primary eclipse, time of periastron passage, rate of apsidal motion, apsidal period, and inclination angle. The 
ARV entries represent offsets between the primary and secondary velocities in each data set (primary minus secondary), and 
an overall difference in th e velocity zero points in the sense CfA minus Fairborn. The physical constants used here are those 
adopted by ITorresI II2010I) . 


5. DIFFERENTIAL PHOTOMETRY AND LIGHT CURVE 
SOLUTION 

An extensive program of CC D photometry was c arried 
out using the NFO WebScope (jGrauer et alJl200^ near 
Silver City, New Mexico, for the purpose of gathering an 
accurate V-band light curve of VSOlMon for analysis. 
A total of 6729 images were obtained over 281 nights 
between 2005 January and 2015 February. Compari¬ 
son stars (TYC 4799-1956-1 and TYC 4799-1098-1, with 
V = 10.31 and 11.21, respectively, both of spectral type 
AO) within 10 arcmin in the same field of view as the vari¬ 
able star were used to determine differential magnitudes. 
The magnitudes were corrected for nightly variations in 
the photometric zero poi nt as we have don e previously in 
similar studies (see, e.g.. lLacv et al.lU008[ l. The variable 
star differential magnitude was referenced to the magni¬ 
tude of the combined light of the two comparison stars 


(variable minus comparisons) in each image. Table[5]lists 
these differential photometric measurements. 

The V-band light curve of this well-detached sys- 
tem was analyzed using t he JKTEBOP applicat i on of 
Southwoi^ (1201111 (see also iPopoer fc Etzellll981l : lEtzell 
198111 . This program can combine the brightness mea¬ 
surements with radial-velocity measurements and even 
historical eclipse timings, but does not currently take 
into account apsidal motion, nor can it use more than 
one radial-velocity set for each component. As mentioned 
earlier, the determination of duj/dt benefits greatly when 
the eccentricity can be constrained by the spectroscopy. 
Consequently, we have elected to apply JKTEBOP to only 
the brightness measurements of V501 Mon, and to treat 
the radial velocities and eclipse timings separately, as 
already described in Section [2] and Section SI The fit¬ 
ted parameters from JKTEBOP and their uncertainties are 
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Fig. 2.— Ephemeris curve for V501Mon from our combined or¬ 
bital fit, along with the eclipse timing measurements. Primary 
measurements are shown with filled circles. 


TABLE 6 

Light-curve solution for V501 Mon. 


Parameter Value 


. 

. 

rA . 

^B. 

k = vb/ta . 

* (deg). 

«A. 

«B. 

yA = yB . 

e. 

UIA (deg). 

is . 

iA . 

. 

isliA . 

Paid (days). 

TMin I (HJD-2,400,000). .. 

(T (mmag). 

N . 


0.738 ± 0.040 
0.15456 ± 0.00057 
0.0838 ± 0.0013 
0.0707 ± 0.0013 
0.843 ± 0.027 
88.022 ± 0.076 
0.60 ± 0.11 
0.51 ± 0.10 
0.25 fixed 
0.1296 ± 0.0054 
231.0 ± 2.0 
0 fixed 

0.648 ± 0.020 
0.352 ± 0.023 
0.544 ± 0.036 
7.02120688 ± 0.00000050 
53401.64906 ± 0.00016 
13.5980 
6729 



TABLE 5 

Differential V-band measurements of 
V501 Mon. 


HJD 

(2,400,000-1-) 

AV 

(mag) 

53376.65788 

2.672 

53376.66208 

2.704 

53376.66422 

2.697 

53376.66633 

2.708 

53376.66839 

2.726 


Note. — This table is available in its en¬ 
tirety in machine-readable and Virtual Observa¬ 
tory (VO) forms in the online journal. A portion 
is shown here for guidance regarding its form and 
content. 

given in Table [ 6 l and are the central surface brightness 
of the secondary ( Jb) in units of the primary, the sum of 
the relative radii (ta +fB)) the radius ratio (fc = tb/ta), 
the inclination angle (i), the linear limb-darkening coeffi¬ 
cients (ma, ub), the eccentricity (e), the longitude of pe- 
riastron of the primary (wa), the sidereal period (Psid); 
and a reference time of primary eclipse (Tmih i). The 
gravity darkening coefficients (?/a, J/b) were assumed to 
be the same for the two components and were held fixed 
at their theoretical values for radiative stars. The mass 
ratio was adopted from our spectroscopic orbit (Table E]). 
The normalized flux contributions of the stars {£a, ^b) 
were computed internally by the program. Experiments 
in which we allowed for third light (£ 3 ) returned values 
not significantly different from zero. The eccentricity and 
oiA from our fit are perfectly consistent with the much 
better determined values in Table E) A graphical repre¬ 
sentation of the fitted model along with the observations 
appears in Figures [SHSl 

6 . ATMOSPHERIC PARAMETERS AND CHEMICAL 
COMPOSITION 

In this section we describe our more detailed study of 
the atmospheric properties of the V501 Mon components, 
including their elemental abundances, that begins with a 
reconstruction of their individual spectra by the method 


- 0.10 



0.10 E_^^_._._._._^_._._._._^^_d 

0.0 0.5 1.0 

Phase 

Fig. 3.— NFO differential V-band observations of V501 Mon, 
shown with our best model fit. Residuals from the fit are shown at 
the bottom. 

of spectral disentangling (iSimon fc SturmI I1994D . The 
method works best with a time series of high-resolution 
spectra more or less ev enly distributed in orbital phase 
(jHensberge et al.l[2n08D . as we have for V501 Mon. In its 
most general form this procedure can solve for the spec¬ 
troscopic orbital elements in addition to the individual 
spectra, so in a sense it is a g eneralization of the tech - 
nique of Doppler tomography (iBagnuolo fc Gied 119911) . 
in which the separation of the component spectra is made 
for a given set of known RVs for each star. In this case 
the orbit is already well known from our fit in Section El 
and the elements were held fixed, i.e., disentangling was 
performed in ‘pure separation’ mode. 

For this work w e used the disentangling code FDBinary 
(|Iliiic et al.ll2004D , which operates in Fourier space based 
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0.10 


-0.04 -0.02 0.00 0.02 0.04 

Phase 

Fig. 4.— Enlargement of Figure[3]around the primary minimum. 



0.42 0.44 0.46 0.48 

Phase 

Fig. 5.— Enlargement of Figure [3] around the secondary mini¬ 
mum. 

on the prescription of iHadrav^ (|1995[ i , with some im¬ 
provements. The Fast Fourier Transform implemented 
in FDB inary gives more flexibility in selecting the edges 
of the spectral segments for disentangling, preserving the 
original spectral resolution. The choice of where to place 
the edges of each segment is a critical step in Fourier- 
type disentangling as these points should be strictly on 
the continuum. Errors incurred in selecting the edges in 
spectra with a high density of lines or with significant 
rotational broadening can add unwanted undulations in 
the disentangled spectra of the components. 

We ran FDBinary on 27 of our TRES spectra of 



Fig. 6 .— Flux ratio for V501 Mon computed from synthetic spec¬ 
tra with parameters near those of the components (Teff = 7500 K, 
logg = 4.0, and [Fe/H] = -1-0.5 for the primary, and Teff = 7000 K, 
log 5 = 4.0, and [Fe/H] = 0.0 for the secondary), and the measured 
radius ratio k = 0.843. The V-band value estimated from our light 
curve fit (^b/^Ai Table is also shown. 

V501 Mon with the highest signal-to-noise ratios (aver¬ 
aging about 22 per pixel), still retaining good phase cov¬ 
erage. Since imperfections in the normalization of the 
echelle spectra are another potential source of spurious 
undulations, we performed the disentangling in relatively 
short spectral segments, usually 100-150 A wide. The 
full spectral range we analyzed is 4000-6760 A, corre¬ 
sponding to about 30,500 pixels. 

In the absence of a significant variation in the relative 
line strengths of the components as a function of orbital 
phase there is an ambiguity in the determination of the 
intrinsic line strengths for each star. In other words, the 
disentangled spectra remain in the common continuum of 
the total light of the binary system, and renormalization 
to their individual continua requires external informa¬ 
tion, such as light ratios inferred directl y from the light- 
curve fit of an eclipsing binary (e .g., iHensberge et al.l 
120001 : iPavlovski fc Hensbergel l2005ll . For VSOlMon we 
have this (Table H]) , but effectively only at a single wave¬ 
length corresponding to the center of the V band. In 
order to estimate the flux ratio at other wavelengths 
we used synthetic spectra based on PHOENIX models 
(jHusser et al.l[2M^ with parameters near those in Ta¬ 
ble [8] below, and we adopted the measured radius ratio k 
from Table[6l Figure |6] shows the result, which is in good 
agreement with the measured light ratio in V. Therefore, 
for renormalization purposes we used a smoothed version 
of this predicted relation. 

Figure [7] shows a sample section of the disentangled 
spectra for the primary and secondary of V50I Mon (sec¬ 
ond and third from the top) on the common continuum of 
the composite spectrum of the binary (top). Also shown 
at the bottom are the same spectra renormalized accord¬ 
ing to the relation described above from Figure [51 The 
signal-to-noise ratios of the reconstructed spectra are 73 
and 39 per pixel for the primary and secondary, respec¬ 
tively, at a mean wavelength of 5500 A. 

Next we used these disentangled spectra to infer the 
stellar parameters and abundances. Effective tempera¬ 
ture and microturbulence were determined from the 
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Fig. 7.— Portions of the disentangled spectra of V501 Mon before 
and after normalization. Prom top to bottom the figure shows one 
of our composite spectra, the reconstructed primary and secondary 
spectra on the common continuum of the binary, and the same two 
spectra after normalization. 

numerous iron lines, and the usual conditions of exci¬ 
tation balance and null correlation of the iron abun¬ 
dance with the reduced equivalent widths, respectively, 
iterating as necessary. Equivalent widths of the Fe I 
and Fe II lines ca refally selected from the line list of 
[Brmitt_e£^lJ_(l201^were measured with the UCLSYN code 
(|Smallev et al.ll20'^ . which was also used for the calcu¬ 
lation of the theoretical spectra. The surface gravities 
of the stars are well known from our spectroscopic and 
light-curve solutions and were held fixed at the values 
reported later in Table IH Initial results clearly indicated 
an enhanced metallicity for the primary consistent with 
earlier indications of an Am nature, with [Fe/H] ~ -1-0.30. 
We therefore adjusted our atmosphere model accord¬ 
ingly, and repeated the analysis. The resulting excita¬ 
tion temperatures and microturbulence velocities based 
on the Fe I lines are Teff = 7540 ± 140 K and = 
3.1 ±0.1 kms“^ for the primary, and = 6960 ±110 K 
and = 1.3 ± 0.1 kms“^ for the secondary. The 
value for the primary is consistent with results for normal 
A stars of this teniperatu re (e.g.. iKunzli fc Nortblll99^ 
iGebran et al.l I2010L1201411 , while that of the secondary 
seems slightly lower than in other stars of its type. 

With both Fe I and Fe II lines being present in the 
spectra of the V501 Mon components, the effective tem¬ 
peratures can in principle also be derived from the con¬ 
dition of ionization balance, in which the strength of the 
singly ionized Fe lines is influenced by log g as well. How¬ 
ever, with logg hxed as we described above, the iron 
abundances from Fe I and Fe II in each star are already 
in excellent agreement (and are A(Fe I) = 7.83 ± 0.07 
and A(Fe II) = 7.80 ± 0.08 for the primary, and A(Fe I) 
= 7.51 ± 0.05 and A(Fe II) = 7.50 ± 0.07 for the sec- 
ondary) 0 , indicating no need to adjust the T^s values 
since they already satisfy ionization balance. 


We use the standard abundance notation in which A(X) = 


Additional estimates of the temperatures were derived 
by fitting the wings of the Hct and H/3 line profiles 
with the STARFiT code (iTamaio et al.ll201]1 : lKolbas et al.l 
l2014fl . which pe rforms an optimizat ion based on a ge¬ 
netic algorithm (jCharbonneaul [199511 . For this we used 
a precalculated grid of theoretical spectra in local ther¬ 
modynamic al equilibrium (for details, see iKolbas et al.l 
12014 l2015fl . The metal lines superimposed on the wings 
were masked out. The surface gravities and projected ro¬ 
tational velocities (see below) were held fixed, and uncer¬ 
tainties in the derived temperatures were calculated with 
a Markov Chain Monte Carlo technique (Pavlovski et ah, 
in prep.). The results from the Ha lines are 7420± 140 K 
for the primary and 6790 ± 180 K for the secondary, and 
from the H/3 lines they are 7480±150 K and 7070±160 K, 
respectively. The latter results are consistent with our 
earlier estimates from the iron lines, while the Ha val¬ 
ues are somewhat lower, as expected from the shallower 
depth of formation of this line. 

Abundances for individual elements were derived using 
lines with the most reliable log 5 / values, and are listed 
in Table [7l The quoted uncertainties take account of the 
intrinsic scatter from the lines used in each case, as well 
as the uncertainties in logg, and ^t- An additional 
contribution has been added to the error budget corre¬ 
sponding to a 1.5% error in the light ratio of V501 Mon. 
T he results in Table [7] are on the solar abundance scale 
of lAsplund et al.l ()2009H . The abundance pattern of the 
two components is compared in Figure [H While the 
secondary is seen to have solar composition within the 
uncertainties, the pattern for the primary is typical of 
Am stars, with a clear underabundance of Ca and Sc, 
an enhancement of the iron group elements, an d a large 
overabundance of Ba and other heavy elements (iPrestonI 

IT^ . 

The Li 6707 A line is clearly visible in the spectrum 
of the secondary star. Line fitting yields an abundance 
estimate of A(Li) = 3.17 ± 0.16, in which the principal 
contribution to the uncertainty is from the error in the 
effective temperature. This abundance is consistent with 
that expected for stars of this mass and temperature (see, 
e.g.. lRamfrez et ani2012[) . On the other hand, we find no 
trace of the Li line in the primary star spectrum. We can 
only place a rough upper limit of about A(Li) = 2.80. 

Finally, we also made an estimate of the macrotur- 
bulence velocities (Crt) and projected rotational veloci¬ 
ties from the disentangled spectra. For the primary we 
obtained Crt = 6.3 ± 1.6 kms“^ and usinz = 18.4 ± 
1.0 kms“^, and for the secondary Crt = 5.3±1.2 kms“^ 
and V sin* = 15.4 ± 0.9 kms“^. 

7. ABSOLUTE DIMENSIONS 

The combination of the spectroscopic and light-curve 
elements yields absolute masses for V501 Mon with rel¬ 
ative errors of only 0.3% and 0.2% for the primary and 
secondary, and radius errors of 1.5% and 1.8%, respec¬ 
tively. Of the temperature determinations in the pre¬ 
ceding section based on the disentangled TRES spec¬ 
tra, we adopt for each star a weighted average of the 
results from the metal lines and from H/3, which are 
= 7510 ± 100 K and = 7000 ± 90 K. These 

log[n(X)/n(H)] + 12, where n{X.) and n(H) are the numbers of 
atoms per unit volume of element X and of hydrogen. 
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TABLE 7 

Abundances from our disentangled TRES spectra of V501 Mon. 




Primary 


Secondary 



A 

X 

Abundance [^/H] 

N 

Abundance [^/H] 

N 

log £0 


11 

Nal 








6.28 

± 

0.05 

+0.04 

± 

0.06 

6 

6.24 

± 

0.04 

12 

Mg I 

7.76 

± 

0.07 

+0.16 

± 

0.08 

3 

7.52 

± 

0.14 

-0.08 

± 

0.15 

4 

7.60 

± 

0.04 

14 

Si I 

7.64 

± 

0.09 

+0.13 

± 

0.09 

11 

7.52 

± 

0.07 

+0.01 

± 

0.08 

5 

7.51 

± 

0.03 

20 

Ca I 








6.33 

± 

0.09 

-0.01 

± 

0.10 

9 

6.34 

± 

0.04 

20 

Ca II 

5.85 

± 

0.08 

-0.49 

± 

0.09 

6 








6.34 

± 

0.04 

21 

Sc II 

2.52 

± 

0.08 

-0.63 

± 

0.09 

3 

3.11 

± 

0.04 

-0.04 

± 

0.06 

6 

3.15 

± 

0.04 

22 

Ti I 

5.40 

± 

0.05 

+0.45 

± 

0.07 

10 

4.93 

± 

0.09 

-0.02 

± 

0.10 

15 

4.95 

± 

0.05 

22 

Ti II 

5.31 

± 

0.06 

+0.36 

± 

0.08 

11 








4.95 

± 

0.05 

24 

CrI 

6.04 

± 

0.09 

+0.40 

± 

0.10 

15 

5.63 

± 

0.06 

-0.01 

± 

0.07 

44 

5.64 

± 

0.04 

24 

Cr II 

6.05 

± 

0.08 

+0.41 

± 

0.09 

9 

5.61 

± 

0.07 

-0.03 

± 

0.08 

19 

5.64 

± 

0.04 

25 

Mn I 

5.60 

± 

0.11 

+0.17 

± 

0.12 

7 

5.36 

± 

0.06 

-0.07 

± 

0.08 

10 

5.43 

± 

0.05 

26 

Fe I 

7.83 

± 

0.07 

+0.33 

± 

0.08 

159 

7.51 

± 

0.05 

+0.01 

± 

0.06 

126 

7.50 

± 

0.04 

26 

Fe II 

7.80 

± 

0.08 

+0.30 

± 

0.09 

18 

7.50 

± 

0.07 

+0.00 

± 

0.08 

14 

7.50 

± 

0.04 

27 

Co I 








5.05 

± 

0.07 

+0.06 

± 

0.10 

12 

4.99 

± 

0.07 

28 

Ni I 

6.69 

± 

0.08 

+0.47 

± 

0.09 

41 

6.22 

± 

0.07 

+0.00 

± 

0.08 

36 

6.22 

± 

0.04 

39 

YII 








2.17 

± 

0.08 

-0.04 

± 

0.09 

6 

2.21 

± 

0.05 

56 

Ba II 

3.94 

± 

0.11 

+1.76 

± 

0.14 

3 








2.18 

± 

0.09 


Note. — Columns list the atomic number, the element and ionization degree, the logarithm of the 
number abundance on the usual scale in which A(H) = 12, the logarithmic abundance relative to the Sun, 
and the number of spectra l lines measured. The last column gives the reference photospheric solar values 
frorn fAsplund et al.l (I20Q9I '). 



jit 

: Primary : 


.. 




Atomic number 

Fig. 8.— Measured abundances for the primary and secondary 
of V501 Mon. The abundance of Ba II for the primary star is not 
displayed, as it is far off the scale (+1.76, atomic number 56). For 
re ference, the shaded ar ea indicates the standard solar composition 
of lAsplund et ahl l|2Q09ll . 

correspond approximately to spectral types of A6 and 
FO. The Ha estimates are not considered for the rea¬ 
son mentioned above. The measured primary/secondary 
temperature difference, ATeff = 510 zb 130 K, is in very 
good agreement with the value 480 zb 140 K inferred di¬ 
rectly from the light curves through the fitted central 
surface brightness para meter J b and the visual abso¬ 
lute flux calibration of IPoppe'il ()1980[) . We assume in 
the following that the chemical composition of the sys¬ 
tem is represented by the measured abundance for the 
secondary, given that the primary is anomalous. Thus, 
[Fe/H] = -1-0.01 ±0.06, based on the numerous Fe I lines. 


The physical properties of the components are col¬ 
lected in Table |51 For the distance estimate we adopted 
an apparent visual magnitude V = 12.320 ± 0.020 
from the AAVSO P hotometric All-sky Survey (APASS; 
IHenden et al.l[ 2 M^ . along with the absolute lum inosi- 
ties and bolometric corrections from lFloweii (|1996f) . The 
interstellar extinction toward VbOlMon, Ay = 0.55 ± 
0 . 12 , was derived by iteratio ns from the r ecent three- 
dimensional reddening map of iGreen et al.1 (1201511 based 
on Pan-STARRS 1 and 2MASS photometry, and as¬ 
suming Ay = 3AE(B — V). The resulting distance, 
D = 1030 ± 80 pc, is very similar to an inde pendent 
calcu lation relying the visual surface fluxes Isee iPoopeil 
1198011 . which gives D = 1050 ± 70 pc. 

Table [5] lists the predicted rotational velocities of the 
stars under the assumption that tidal forces have syn¬ 
chronized them with the orbital motion at periastron 
(uperisinj), and the alternate assumption that they are 
pseudo-synchronized (upsync sin i : lmhl[TMM 1 Slight dis¬ 
agreements between our three estimates of v sin i from 
Sections [3] and [5] prevent us from distinguishing between 
the two possibilities for V501 Mon. In particular, the two 
non-independent determinations from our CfA/TRES 
spectra differ by about twice their combined errors for 
the secondary component. Weight-averaging them, and 
then combining the result with the Fairborn estimates 
for each star yields mean vsini values of 17.1 ± 0.7 and 
12.8 ± 1.2 kms“^, although we suspect systematic errors 
are not negligible. 

8. COMPARISON WITH STELLAR EVOLUTION MODELS 

Figure [3] presents a simultaneous comparison of the 
masses, radii, and temperatures of the V501 Mon com¬ 
ponents wit h stellar evolu t ion models from the Y onsei- 
Yale series (|Yi et al.ll200H [Pemaraue et al.l[2003l . The 
solid lines correspond to evolutionary tracks for the ex¬ 
act masses we measure, beginning at the zero-age main 

^ Both scenarios also assume the stars’ spin axes are parallel to 
the orbital axis. 
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TABLE 8 

Physical properties of V501 Mon. 


Parameter Primary Secondary 


Mass (Mq). 

. 1.6455 

± 

0.0043 


1.4588 

± 

0.0025 

Radius (R©). 

1.888 

± 

0.029 


1.592 

± 

0.028 

logg (cgs). 

4.103 

± 

0.013 


4.199 

± 

0.016 

Temperature (K). 

7510 

± 

100 


7000 

± 

90 

logL/L©. 

1.007 

± 

0.027 


0.743 

± 

0.043 

. 

0.03 

± 

0.10 


0.03 

± 

0.10 

Mbol (mag)'^. 

2.215 

± 

0.067 


2.87 

± 

0.11 

Mv (mag). 

2.18 

± 

0.12 


2.84 

± 

0.15 

Fv^ . 

. 3.8749 

± 

0.0057 


3.846 

± 

0.010 

Av (mag). 



0.55 

± 

0.12 



Distance (pc)*^. 



1030 

± 

80 



m — M (mag). 



10.06 

± 

0.16 



Uperi sin 2 . 

18.0 

± 

0.3 


15.1 

± 

0.3 

'i^psync sin i . 

15.1 

± 

0.2 


12.7 

± 

0.2 

Measured v sin i ®. 

16.6 

± 

1.0 


11.2 

± 

2.0 

Measured v sin. 

18.4 

± 

1.0 


15.4 

± 

0.9 

Measured usini®. 

16.5 

± 

1.0 


12.4 

± 

1.0 

(kms-l). 

3.1 

± 

0.1 


1.3 

± 

0.1 

Crt (kms“l) . 

6.3 

± 

1.6 


5.3 

± 

1.2 

[Fe/H] (dex). 



+0.01 

± 

0.06 



A(Li 6707 A) (dex)... 

< 

2.80 


3.17 

± 

0.16 


® Bolometric corrections from IFIoweii Ill99fili . 

^ Uses = 4.732 fsee lTorre3l20ini . 

^ Visual absolute flux IIPopDei1ll98^ . 

^ Relies on the luminosities and bolometric corrections. 

® Based on the original (composite) TRES spectra, determined 
from cross-correlation grids. 

^ Based on the disentangled TRES spectra, 
s Based on the Fairborn spectra. 

sequence (ZAMS). The (solar-scaled) chemical compo¬ 
sition of the models has been adjusted to reproduce 
the observations, and the best match is achieved for 
[Fe/H] = —0.07, which is not far from the measured 
value of [Fe/H] = -1-0.01 ± 0.06. An isochrone for this 
metallicity and the best-fit age of 1.13 Gyr is shown with 
a dashed line. The stars are seen to be slightly evolved 
from the ZAMS. The radii and temperatures are shown 
separately as a function of mass in Figure [TUI and com¬ 
pared against the same best-fit model isochrone. The 
agreement is very good. 

The difference between the best-fit metallicity from the 
models and the measured composition of V501 Mon cor¬ 
responds to about 1.3 times the uncertainty. If due to 
observational error in the temperatures and/or compo¬ 
sition, one possible source for the small discrepancy is a 
bias in the light ratio we have adopted for the renormal¬ 
ization of the disentangled TRES spectra. The adopted 
dependence of the light ratio with wavele ngth relies on 
synthe tic spectra for normal stars from iHusser et al.l 
(|2013f) . which for the primary may differ from the true 
spectrum in subtle ways because of the Am nature of 
that star. 

9. DISCUSSION AND CONCLUDING REMARKS 

The high-precision masses and radii derived here for 
detached, double-lined eclipsing binary VSOlMon place 
it in the elite gr oup of systerns wit h the best determined 
properties fe.g.. [Torres et al.ll201flh . Our detailed chemi¬ 
cal analysis adds special value as fewer than three dozen 
of the systems with masses and radii known to better 
than 3% have been subjected to this type of study. The 
discovery that the primary is a metallic-line A star makes 
it an important object that may help us understand more 



Fig. 9. — Stellar evolution models from the Yonsei-Yale series 
compared against the observations for V501 Mon. Evolutionary 
tracks (solid lines) are shown for the exact masses we measure, and 
for the metallicity that produces the best match, which is near the 
measured value. The very small mass errors mean that the uncer¬ 
tainty in the placement of the tracks is very small (indicated in the 
plot with a barely visible shaded region around each track). The 
best-fit isochrone shown by the dashed line has an age of 1.13 Gyr. 




Fig. 10.— Mass-radius and mass-temperature diagrams show¬ 
ing the reference isochrone from Figure 0 ([Fe/H] = —0.07, 
1.13 Gyr). The error boxes represent the observational uncertain¬ 
ties for V501 Mon. 
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about this phenomenon. So far as we are aware, only four 
other eclipsing binary systems with similarly well mea¬ 
sured masses and radii have been published in which one 
or both components are Am/Fm stars, and which also 
have d etailed element-by-elein ent abundance analyses: 
dAur (iLvubimkov et al.(il99fifl . S W CMa and HW CM a 
(|Torres et al.l 2012I1 . and YZ Cas (iPavlovski et al.l[2014ll . 

Figure im shows all well-studied (ctm and less than 
3%) A- and F-type stars in the field that are in the 
te mperature range of V501 Mon from the compilation 
of [Torres et al.l ( 20101 1 , with additions from the recent 
literature. Metallic-line stars are indicated with filled 
circles, and VSOlMon is represented with squares. Nu¬ 
merous studies in the literature have investigated the 
properties of Am/Fm stars in the field or in clusters. 
A connection has been sought between abundance pat¬ 
terns in these objects and various characteristics such 
as binary orbital period or ec centricity, rotation , tem¬ 


perature, age, and others (e.g.. iBudailHQQGL 11997 


kunzli fc Northl 119981 lAM l2000l: Fehovcfk et al 


1999; 


2004; 


North fc Debernardil 120041 iPrieur et al.ll2006L and refer¬ 
ences therein), but generally without an accurate knowl¬ 
edge of the basic stellar properties (mass, radius) or the 
age, which may be important if the anomalies depend 
sensitively on them. Only seven of the points in Fig¬ 
ure [TT] (from VSOlMon and the four systems mentioned 
earlier) have detailed photospheric abundances and many 
more may be needed before firm correlations can be¬ 
gin to be discerned as a function of mass, radius, or 
perhaps age. For now, the figure shows that the cur¬ 
rently well-measured Am/Fm systems in the field are 
confined to the age range from about 0.4 Gyr to 1.5 Gyr 
(indicated with dashed lines), whereas normal A or F 
stars can be younger or older. The upper boundary 
also correspo nds to a roughly const ant temperature of 
6700-6800 K. iKunzli fc Northl (|1998[1 found similar lim¬ 
its for field stars. Younger Am stars do of course exist 
in open c lusters, such as the handful of examples in the 
Pleiades (lAbt fc LevatolfTOTSl iBurkhart fc GouDrvlll997l : 


iHui-Bon-Hoa fc Alecianlll99^ iGebran fc Monied 1200811 

and perhap s even among pre-main-sequence s tars such 
as AKSco (jHerbigI [19601: lAndersen et al.lll989l l with an 
age of only ^18 Mvr (iGzekala et al.ll2015ll . 

We detect Li in the spectrum of the secondary of 
VSOlMon but not in the Am primary. Some studies 
have suggested that the Li abundance in Am stars is 
lower than in normal A sta rs by about a factor of 3 (e.g., 
IBurkhart fc Couprvl I2000l) , while others have found no 


significant difference fe.g.. lCatanzaro fc Rip enfc 01411 . In 

the 1.1 Gyr system V501 Mon the equal-age components 
do differ by at least a factor of 2, though it is unclear 
whether this may be more related to the temperature 
difference (ATeg ~ 500 K in this case). 

Finally, the measurement of apsidal motion in eclips¬ 
ing binary systems has long been used to test models 
of the internal structure of stars (for a recent review, 
see iClaret fc Gimenedl201^ . specifically, the degree of 
mass concentration towards the center. In addition to 
the classical contributions to dui/dt from tidal and rota¬ 
tional distortions, i n a few systems th e contribution from 
General Relativity (|Levi-Givitalll937ll is very important. 
This happens to be the case for V501 Mon, for which we 
calculate (da;/dt)GR = 0.00032 deg cycle“^. Although 


the total measured apsidal motion is still very uncer- 



Fig. 11.— Field eclipsing binaries (circles) with well deter¬ 
mined masses and radii in the temperature range of V501 Mon 
(squares). Solid lines connect components of the same binary, 
and metallic-line stars are indicated with filled symbols. Solar- 
metallicity isochrones from the Yonsei-Yale series (0.1-13 Gyr) are 
shown with dotted lines for reference. The dashed lines delimit 
the age range 0.4-1.5 Gyr within which the Am/Fm stars in this 
sample fall. 


tain [doj/dt = 0.00045 ± 0.00024 deg cycle~^) , the GR 
term appears to dominate, as anticipated by iGimene^ 
(jl985fl . nominally contributing about 70% to the total 
amount. The internal structure constants for the pri- 
mary and se condary of V501 Mon from the models by 
I Claret] (|2004il . logfc 2 = —2.51 and —2.42, lead to a pre¬ 
dicted total apsidal motion rate (including the GR term) 
of {duj/dt)tot = 0.00046 deg cycle“^ that is very close to 
the measured value. 
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